Microscopic calculation of proton capture 

reactions in mass 60-80 region and its 

astrophysical implications 

Chirashree Lahiri and G. Gangopadhyay^ 
Csl '. Department of Physics, University of Calcutta 

92 A.P.C. Road, Kolkata 700009, India 
^e-mail : ggphy@caluniv.ac.in 

July 25, 2011 

-J , Abstract 

5 . Microscopic optical potentials obtained by folding the DDM3Y 

interaction with the densities from Relativistic Mean Field approach 
have been utilized to evaluate S-factors of low-energy (p, 7) reactions in 
^ . mass 60-80 region and to compare with experiments. The Lagrangian 

Q\ I density FSU Gold has been employed. Astrophysical rates for impor- 

Q ■ taut proton capture reactions have been calculated to study the be- 

,—1- . haviour of rapid proton nucleosynthesis for waiting point nuclei with 

, I [ mass less than A = 80. 
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1 Introduction 

Relativistic Mean Field (RMF) calculations have proved to be very 
5^ ■ successful in describing different features of nuclei. This method has 
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been used to study binding energy of the ground state and vari- 
ous excited states, deformation, charge radii, density profile and nu- 
clear halo, etc[rj. Particularly, the success in reproducing the density 
has motivated various microscopic calculations that generate nucleon- 
nucleus and nucleus-nucleus potentials for study of proton and alpha 
radioactivity and scattering (See [21 [3l S] and Refs. therein). In the 
present work, we use RMF densities to produce microscopic optical 
potentials [5] to study proton capture reactions in mass 60-80 region. 
Proton capture reactions at very low energy play a very important 
role in nucleosynthesis. Particularly rapid proton capture (rp) process 



in explosive nucleosynthesis is a basic ingredient in driving the abun- 
dance along the N = Z line[6l [3- As this process has to overcome a 
large Coulomb barrier it can occur only at a higher temperature range. 
For example, X-ray bursts provide a large flux of protons at peak tem- 
peratures around 1-2 GK and are expected to play a significant role 
in the creation of nuclei up to mass 110. 

The rp-process proceeds along the N = Z line in mass 60-80 region. 
In nature, the important proton capture reactions usually involve cer- 
tain nuclei as targets which are not available to us. Hence, experi- 
mental information is difficult, if not impossible, to obtain, at least in 
near future. In such a situation, one has to rely on theory for the reac- 
tion rates. Rauscher et al. have extensively calculated reaction rates 
and cross sections in a global approach [51 [9]. They have commented 
that statistical model calculations may be improved by using locally 
tuned parametrization of nuclear properties such as optical poten- 
tial. However, they prefer a global approach to predict astrophysical 
rates for experimentally inaccessible nuclei. Cross sections have been 
calculated for proton capture reactions in mass 60-80 regions using 
a semimicroscopic optical potential in the local density approxima- 
tion using phenomenological density prescriptions. However, far from 
the stability valley, these prescriptions may not represent the actual 
densities very well leading to considerable uncertainty in the reaction 
rates. Very often, the reactions rates are varied by a large factor to 
study their effect. For example, SchatzflO] varied the rates of certain 
reactions by a factor of hundred. Obviously, this makes the results 
uncertain to some extent. 

A fully microscopic calculation may be used to estimate the rates 
to reduce the above uncertainty. A consistent framework for calcula- 
tion may be constructed based on microscopic densities. This has the 
advantage of extending it to unknown mass regions. In the present 
work, we have tried to calculate the reaction rates from a purely mi- 
croscopic model, i.e. RMF. We have already mentioned that RMF is 
particularly suitable to describe nuclei far away from the stability val- 
ley where experimental information is scarce. A microscopic optical 
potential obtained by folding an appropriate microscopic NN inter- 
action is expected to be more accurate and may do away with the 
necessity of any arbitrary variation in the reaction rates. 

However, even in microscopic optical model, there often remain 
certain parameters which can be fixed only after comparison with 
experiment. We have compared the results for a number of reactions 
in the mass region ^ = 60 — 80 for which experimental informations 
are available. This has helped us in determining a set of parameters 
for this mass region. 



Once the parameters have been fixed, we employ them to calculate 
the rates of a number of astrophysically important proton capture re- 
actions. Certain N = Z nuclei having the highest abundance in an 
equilibrium in a chain are called the waiting points^ for the chain. 
These nuclei have negative or small positive Q-values for proton cap- 
ture. An equilibrium between the (p, 7)/(7,p) processes is established 
and the rp process may have to wait for beta-decay or a-capture to 
proceed to heavier nuclei. Certain N = Z waiting point nuclei with 
A < 80, viz. ^^Zn, ^^Se, '''^Kr, and ^^Sr have long half lives, their total 
lifetime being large compared to the time scale of typical X-ray bursts 
(10-100 sec). Thus, they may produce a bottleneck in the rp-process 
that would slow down the rate of hydrogen burning and necessitate 
extended burst tails unless two proton capture can reduce these half 
lives and bridge the waiting points. X-ray burst model calculations 
are therefore particularly sensitive to the rates of proton capture for 
these nuclei. We have used the microscopic approach, used in the 
present work, to calculate the rates with an aim to study the bridging 
of the waiting point nuclei. 

2 Method 

As already mentioned, experimental (p, 7) rates in many nuclei in- 
volved in the rp-process are not available as they are unstable. Hence, 
theory remains our sole guide. A microscopic optical model calcula- 
tion based on theoretical mean field densities are expected to provide 
us with reliable rates. With this in mind, we have calculated the nu- 
clear density profiles in a RMF approach using the Lagrangian density 
FSU Gold [11]. This density contains two additional non- linear meson- 
meson interaction terms, whose main virtue is a softening of both the 
Equation of State of symmetric matter and symmetry energy. As a 
result, the new parametrization becomes more effective in reproducing 
quite a few nuclear collective modes, namely the breathing modes in 
^^Zr and ^'^^Pb, and the isovector giant dipole resonance in ^'^^PbpTj. 
We have found this Lagrangian density to be very useful in explaining 
various other properties that depend on the nuclear density, such as 
nuclear decay and reaction [21 [3lll]. 

As the important quantity in calculating proton capture cross- 
section is the nuclear density as a function of radius, the calculations 
have been carried out in co-ordinate space assuming spherical symme- 
try. Pairing has been introduced under the BCS approximation using 
a zero range pairing force of strength 300 MeV-fm for both protons 
and neutrons. The RMF+BCS equations are solved under the usual 
assumptions of classical meson fields, time reversal symmetry, no-sea 



contribution, etc. The details of the calculation may be obtained from 
Bhattacharya et al. [12l [13] . 

Microscopic optical model potential is usually obtained by folding 
an effective interaction, derived from the nuclear matter calculation, in 
the local density approximation, i.e. by substituting the nuclear mat- 
ter density with the density distribution of the finite nucleus. In the 
present work, the density dependent M3Y (DDM3Y) effective inter- 
action [m [151 [121 [TTj has been utilized for this purpose. This was 
obtained from a finite range energy independent G-matrix elements of 
the Reid potential by adding a zero range energy dependent pseudopo- 
tential and introducing a density dependent factor. The interaction is 
given by 

v{r,p,E)=t'^'^''{r,E)g{p) (1) 

where E is incident energy and />, the nuclear density. The i^^^v 
interaction is given by 

tM3Y ^ 7999 e _ 2134^-— + JooiE)6{r) (2) 

4r 2. or 

where Joo{E) is the zero range pseudo potential, 

Joo{E) = -276 ( 1 - 0.005^ j MeV/m^ (3) 

and g{p) the density dependent factor, 

g{p) = C{l-bp'/^) (4) 

The constants in the last equation have been obtained from nuclear 
matter calculation |17j as C = 2.07 and b = 1.624 fm'^. We have used 
this form in our calculation keeping the above parameters unchanged. 
Since nuclear matter-nucleon potential does not include a spin- 
orbit term, the spin-orbit potential from the Scheerbaum prescription [T8 
coupled with the phenomenological complex potential depths X^so and 
Xwso has been used. 

K^i^) = (^-"so + ^>'wso)--^{^Pp{n) + -^Pn{p)) (5) 

The depths are functions of energy, given by 

Kso = 130exp(-0.013£;) + 40 

X^so = -0-2{E - 20) 

where E is in MeV. These standard values have been used unaltered 
in the present work. 



Cross-section and astrophysical rates are calculated in the Hauser- 
Feshbach formalism using the computer package TALYS1.2[19j. The 
present method has already been applied in mass 60 region in Ref.[3]. 
In the present approach, though, a slightly different normalization has 
been applied as described later. 

For the calculation of proton capture at waiting points, a small net- 
work has been designed which includes the following processes. The 
waiting point nucleus with Z = N, which acts as a seed, may capture 
a proton. The resulting nucleus, with Z = N + 1, may either capture 
another proton or undergo photodisintegration emitting a proton to 
go back to the seed nucleus. The nucleus with Z = N-\-2 may also un- 
dergo photodisintegration. In addition, all the three nuclei mentioned 
above may undergo /3-decay. The photodisintegration rates at differ- 
ent temperatures have been calculated from the proton capture rates 
using the principle of detailed balance. The density has been taken 
as 10^ gm/cm'^ unless otherwise mentioned. The proton fraction has 
been assumed to be 0.7. 

One of the difficulties in the calculation is the unavailability of 
sufficiently accurate experimental Q-values in most cases. The pho- 
todisintegration rate is exponentially dependent on the Q-value. In 
some of the Z = N + 1 nuclei, experimental binding energy values are 
either not available, or have very large errors. In absence of experi- 
mental values, we have used the Q-values and their errors adopted in 
Refs. |10] and [20j for the network calculation. TALYS code uses the 
Dufio-Zuker|2H [22] formula for masses for which experimental values 
are not available. These have been kept unchanged as the small dif- 
ference between the values predicted by the formula and the values 
adopted in pH] and [2D] do not significantly affect the reaction rates. 
The Q-values used for the calculation are indicated in table 1. 

The measured half life values for beta-decay have been taken from 
the compilation by Audi et al. [26] except in the case of the case of 
^^As. For this nucleus, a measured value of 0.128 second from Ref. 
|27] has been assumed. If measured values are not available, we have 
adopted the values from the work by Moller et al. [28] . 

3 Results 

This section is subdivided in three parts. In the first part, we dis- 
cuss the results of RMF calculations very briefly. We concentrate on 
the results on the nuclei where experimental measurements on (p, 7) 
reactions have been performed. In the second part, we compare the 
available reaction data with the microscopic calculations obtained in 
the procedure described above. In the third subsection, we discuss the 



proton capture reactions at the waiting points. 

3.1 Relativistic mean field calculations 

We present a very brief summary of the resuhs for only those nu- 
clei for which experimental measurements of {p, 7) reaction are avail- 
able. In table 2 we compare the theoretical results for binding energy 
and charge radius with experimental information wherever available. 
Charge radius has been chosen for comparison as it provides the sim- 
plest measure of the charge density distribution in the nucleus. Some 
of the binding energy results were presented in Ref. 0] earlier. The 
binding energy values from the mean field approach have been cor- 
rected using the formalism developed in Ref. [291 SO]- The experi- 
mental binding energy and charge radius values are from Refs. |2L)| 
and [32], respectively. 

Charge radius has been calculated from the charge density of the 
nucleus. The charge density, in turn, has been obtained from the point 
proton density pp by taking into account the finite size of the proton. 
The point proton density is convoluted with a Gaussian form factor 
5(r), 

Pchir) = I epp{r')g{r - r')dr' (6) 

5(r) = (aV^)-3exp(-rVa2) (7) 

with a = 0.8 fm. We plot in Fig. [T]the calculated charge density for 
^^Ni and ^^Zn as representatives of our results. Experimental mea- 
surements from Wohlfahrt et al.[61\ are shown as filled circles. One 
can see that the theoretical and experimental values agree very well, 
particularly at larger radii values, which is the region expected to con- 
tribute to the optical potential at low projectile energy. Other nuclei 
also show similar agreement. 

One can see that the charge radius values are also reproduced to 
considerable accuracy. There were two obvious typographical errors 
in the theoretical charge radii values of ^^'^"^Ni in Ref. [3]. One can 
see that the experimental quantities are adequately described in the 
present formalism. 

3.2 S-factors of (j?, 7) reactions 

Cross-sections of low-energy {p, 7) reactions, for which experimental 
data are available in the A = 60—80 region, have been calculated. The 
energy relevant to the rp-process in this mass region lies between 1.1 



to 3.6 MeV. As the cross-section varies very rapidly at such low energy, 
a comparison between theory and experiment is rather difficult. The 
usual practice in low energy nuclear reaction is to compare another 
key observable, viz. S-factor. It is given by 

S{E) = Ecj{E)e^^'^ (8) 

where E is the energy in centre of mass frame in KeV, (j{E) indicates 
reaction cross-section in barn and r/ is the Sommerfeld parameter with 

27rr? = 31.29ZpZm/-^ (9) 

V E 

Here, Zp and Zt are the charge numbers of the projectile and the 
target, respectively and // is the reduced mass (in amu). The quan- 
tity S-factor varies much more slowly than reaction cross-sections as 
the exponential energy dependence of cross-section is not present in 
it. For this reason, we calculate this quantity and compare it with 
experimentally extracted values. 

As already pointed out, our calculations, being more microscopic, 
are more restricting. Yet, the rate depends on the models of the level 
density and the El gamma strength function adopted in the calcula- 
tion of cross sections. Phenomenological models are usually fine tuned 
for nuclei near the stability valley. Microscopic prescriptions, on the 
other hand, can be extended to the drip lines, which is a requirement 
for calculating the rates at the waiting points. Hence, microscopic ap- 
proach has been assumed for all nuclei. We have calculated our results 
with microscopic level densities in Hartree-Fock (HE) and Hartree- 
Fock-Bogoliubov (HFB) methods, calculated in TALYS by Goriley 
and Hilaire, respectively. For El gamma strength functions, results 
derived from HF-I-BCS and HFB calculations have been employed. 
All these options are available in the TALYS data base. 

The real part of the potential has been obtained by normalizing 
the folded DDM3Y potential by a factor of 0.7, while the imaginary 
part, by a factor of 0.1, so as to explain the S-factors obtained in the 
above experiments. In Ref. [3], we used a slightly different normaliza- 
tion. There, for a good fit, we also used another parameter, Gnorm, for 
normalizing the gamma-strengths. However, we found that this pa- 
rameter varies from nucleus to nucleus. If the present calculation has 
to be extended to unknown nuclei, the approach is clearly inadequate. 
To overcome this difficulty, we have used a different normalization in 
the present work as stated above. 

In Figs. 2-4, we plot the results of our calculations against S-factors 
extracted from experimental cross section values. The target nuclei are 
indicated in the figures. In ^^'^^Ni, the experimental values are from 



Tingwell et a/. [33] and Sevior et a/.[3l], respectively. For ^^'^^Cu, the 
experimental results are from Ref. [34j (empty circles) and the Ph.D. 
thesis of Qiang[35] (filled circles). In Fig. [3l the experimental values 
for ^^Zn are extracted from Refs.[36j (empty circles), [37] (filled circles) 
and [38] (empty triangles). The results for ^^Zn are from Skakun et 
al. [3H] while those from ^^'^'^Zn are from Refs. \6<6\ I39j. respectively. 
The results for Ge are taken from Kiss et al.[l0]. S-factors for ' Se 
are from Gyurky et aL[lTj and Krivonosov et a/.[36j, respectively. S- 
factors for ^*^Se(p,7) reactions are also from Refs. [36] (empty circles) 
and [H] (filled circles). In some instances, the actual numerical values 
have been obtained from the website of the National Nuclear Data 
Center [52] of the Brookhaven National Laboratory. 

We find that, except in the case of Se the results for HF plus 
BCS and HFB calculations describe the S-factors reasonably well. 
They both reproduce the general trend and, except for a few isolated 
points, come very close to experimental values. In many of the above 
measurements, the errors were not given. Some of the measurements 
are also old. Taking all these facts in to consideration, it is easy to 
see that the HFB method scores particularly well. Interestingly, more 
recent experimental data, whenever available, show better agreements 
with HFB results than previous measurements. The HFB approach 
also is known to work well for nuclei away from the stability valley. 
So, for later calculations, we have employed the level density and El 
gamma strength values from HFB calculation. The normalization of 
the potential, as described above, has also been assumed. 

3.3 Waiting point nuclei in rp process 

Effective half life values for seed nuclei ^^Ge, ®^Se, ^^Kr, and ^^Sr are 
shown in Figs. 5-8 using continuous lines. For comparison, we have 
also plotted the results calculated from the rates in Rauscherei a/.^ 
by dash-dotted lines. The effects of the uncertainties on the half life 
values in the Q-values have been indicated in the figures by dotted 
lines. 

One can see that the effective half life of the waiting point nuclei 
indeed gets reduced by two proton capture, reaching the minimum 
around 1-2 GK temperature. However, the rate of change is different 
at different waiting points. We need to check whether the reduction 
is sufficient for nucleosynthesis to bridge the waiting points and to 
proceed along the Z = N line. 

Fig. [5] shows the change in the effective half life of ^^Ge in explo- 
sive hydrogen rich environment. The half life decreases and possibly 
goes to a value substantially less than ten seconds, the minimum du- 



ration of an X-ray burst. However, one sees that the uncertainty in 
mass measurement prevents one from reaching any firm conclusion. 
Depending on the actual value of the masses, it may even be possible 
that a burst of the order of ten seconds cannot bridge this waiting 
point effectively. 

In Figs. and [7] one can see that effective half lives of ^^Se and 
^^Kr are not affected significantly by rapid proton processes. This 
conclusion remains unchanged even if the errors in Q- values are taken 
into account. Thus, the most likely scenario is that the rp-process 
path shifts to more stable nuclei and proceed along a different path. 
Possibilities of bridging these two points in astrophysical environment 
have been briefly discussed later in this work. In the case of ^^Sr, 
on the other hand, the half life decreases by more than a factor of 
two. It is clear that this particular waiting point is easily bridged by 
rp-process. 

Is it possible for the waiting points at ^*Se and ^^Kr to be bridged 
in environmental conditions other than assumed in the calculation? If 
the two proton capture rate exceeds that of beta-decay at the waiting 
point, we may consider it to be effectively bridged by rapid proton 
capture. In Fig. 9 we show the densities above which the waiting 
points at ^^Ge, ^^Se, ^^Kr, and ^^Se are bridged, as a function of tem- 
perature. We see that for a density of lO^gm/cm^, the first of the 
above waiting points are bridged over a large temperature range. The 
waiting point at ^^Sr has a narrower temperature window but can def- 
initely be bridged at the density. The point at ^^Se requires slightly 
higher density which may perhaps be available in astronomical envi- 
ronments. However, the waiting point at ^^Kr presents a completely 
different picture. One sees that a density in excess of 10^ gm/cm^ is 
required for effectively bridging this waiting point. Such a high den- 
sity is not expected even in X-Ray burster environment. Thus this 
waiting point is expected to stall the rp-process nucleosynthesis and 
shift it towards the stability valley. 

We note that the rates from Rauscherei al.[Q\ produce nearly iden- 
tical results except in the case of ^^Kr. The rates used in the present 
calculation are slightly smaller, a fact that gives rise to a slight de- 
crease in the temperature where the half life has a minimum. As our 
calculation is microscopic in nature, while at the same time depend- 
ing on local density approximation, we expect the results obtained in 
the present work to be more reliable. In all the cases studied above, 
including that at Kr, the conclusions arrived at the present work do 
not change to any appreciable extent in case the rates from Ref. [5] 
are adopted. This is to be expected as both the present work and Ref. 
[8] use the Hauser-Feshbach formalism. 



4 Summary 



In summary, RMF calculations has been performed in various nu- 
clei in mass 60-80 region using FSU Gold Lagrangian density. Cross 
sections for (p, 7) reactions in nuclei with 60 < ^ < 80 have been cal- 
culated using a microscopic optical potential obtained by folding the 
DDM3Y effective interaction with the theoretical nuclear densities. 
The parameters employed in the microscopic optical potential have 
been fixed after comparison of experimental S-factors with calculated 
ones. We have found that the level density and El gamma strength 
obtained form HFB calculation are eminently suited to describe the 
observed S-factors. Astrophysical rates for proton capture and photo- 
disintegration rates have been calculated in rp-process waiting point 
nuclei with 60 < A < 80 using the same approach. The possibility of 
bridging the waiting points by two proton capture has been investi- 
gated. We see that unlike the other waiting points studied, the one at 
^^Kr is unlikely to be bridged by two proton capture. 
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Table 1: Q- values adopted in the calculation 
Reaction Q-values(MeV) 

^^Ge{p,-ff^As -0.255± 0.104 [23j 

6^As(p, 7)6 '^Se 2.350 ± 0.200 [23j 

'58Se(p,7)69Br -0.785 i°:°34 [24j 

''^Br(j9,7)™Kr 2.605± 0.16 [ID] 

^2Kr(p, 7)^3j^b -0.71 ± 0.10 [25] 

^3Rb(V,7)^^Sr 2.18 ± 0.14 [25] 

^6Sr(p,7)^^Y -0.050 ± 0.072 [20] 

^^Y(p, -fY^Zr 2.087 ± 0.507 (20] 



Table 2: B.E. and radius values in selected nuclei 



Nucleus 


Binding 


Energy (MeV) 


Charge 


radius (/m) 




Theo. 


Exp. 


Theo. 


Exp. 


62Ni 


544.71 


545.26 


3.854 


3.841 


64Ni 


561.97 


561.76 


3.868 


3.859 


63Cu 


551.17 


551.38 


3.889 


3.883 


65CU 


569.43 


569.21 


3.908 


3.902 


e^Zn 


558.815 


559.098 


3.923 


3.929 
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Figure 1: Calculated charge density in ^^Ni and ^^Zn (solid lines) compared 
with experimental measurements (filled circles) from Ref. 
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Figure 2: S-factors extracted from experimental measurements compared 
with theory for ^"'^^Ni and ^^'^^Cu. Sohd and dashed hues indicate respec- 
tively the results of the HF+BCS and HFB approaches for level density and 
El gamma strength. Ecm is centre of mass frame energy. 
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Figure 3: S-factors extracted from experimental measurements compared 
with theory for 64,66,67,682^_ ggg caption of Fig. |2]for details. 
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Figure 4: S-factors extracted from experimental measurements compared 
with theory for ™Ge and ^^''^^''^^Se. See caption of Fig. |2]for details. 



17 



o 

i-H 

X 

o 



10^ 



10^ 



10 



10' 










0.8 



1.0 1.2 1.4 1.6 
Temperature (GK) 



1.8 



Figure 5: Effective half life values of ^^Ge as a function of temperature. The 
solid line represents the results of our calculation while the dashed lines mark 
the two extremes for the errors in the Q- values of the reactions involved. The 
dash dotted line shows the results obtained using the rates from [8]. 
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Figure 6: Effective half life values of ®^Se as a function of temperature. See 
caption of Fig. [5] for details. 
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Figure 7: Effective half life values of ^^Kr as a function of temperature. See 
caption of Fig. [5] for details. 
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Figure 8: Effective half life values of ™Sr as a function of temperature. See 
caption of Fig. [5] for details. 
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Figure 9: Densities above which the various waiting points are found to be 
to be effectively bridged at different temperatures. The waiting points are 
^■^Ge (continuous curve), ^^Se (dotted curve), '^^Kr (dash-dotted curve) and 
^^Sr (dashed curve). 
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